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ABSTRACT 

A key argument in favor of orientation based unification schemes is the finding that 
among the most powerful 3CRR radio sources the (apparent) median linear size of 
quasars is smaller than that of radio galaxies, which supports the idea that quasars 
are a subset of radio galaxies, distinguished by being viewed at smaller angles to the line 
of sight. Recent measurements of radio sizes for a few other low frequency samples are, 
however, not in accord with this trend, leading to the claim that orientation may not 
be the main difference between radio galaxies and quasars. We point out that this "in- 
consistency" can be removed by making allowance for the temporal evolution of sources 
in both size and luminosity, as inferred from independent observations. This approach 
can also readily explain the other claimed "major discrepancy" with the unified scheme, 
namely, the difference between the radio luminosity-size correlations for quasars and 
radio galaxies. 

Subject headings: galaxies: active — galaxies: jets — galaxies: quasars: general — radio 
continuum: galaxies 
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1. Introduction 

According to a widely discussed unified scheme for 
powerful F-R II extragalactic radio sources (with lu- 
minosities above ~ 10 25 ' 5 W Hz^ 1 at 1 GHz, where 
we take Ho = 50 km s _1 Mpc _1 and qo — 0.5 through- 
out this Letter), narrow-line radio galaxies (RGs) are 
identified as quasars (QSRs), or even blazars, when- 
ever their principal axis happens to be oriented within 
a certain critical angle, ip, from the line-of-sight (for 
recent reviews, see Antonucci 1993; Urry & Padovani 
1995; Gopal-Krishna 1995, 1996). In this model, 
the parsec-scale nuclear core of such sources, consist- 
ing of a compact central engine ejecting relativistic 
jets of non-thermal continuum emission, and a broad 
line region (BLR), is believed to be surrounded by 
a dusty torus. In the case of RGs, the torus is be- 
lieved to obscure the core in the visible through soft 
X-ray bands, so that the BLR is not directly visi- 
ble. Strong evidence for such tori comes from the de- 
tection of the BLR in the (scattered) polarized light 
(e.g., Antonucci & Miller 1985; Cimatti et al. 1993; 
Draper, Scarrott & Tadhunter 1993). Various evi- 
dences for the relativistic beaming aspect of this pic- 
ture include: extremely rapid variability of blazars 
over all bands; apparent superluminal motion in the 
bright radio cores; larger ratios of core to total ra- 
dio emission, / c , in QSRs than in RGs; correlation 
of f c with the polarized optical nuclear continuum 
and its anti-correlations with apparent radio linear 
size, the symmetry of core-to-lobe separations, and 
the equivalent width of the [O II]A3727 line (Urry & 
Padovani 1995; Gopal-Krishna 1995, and references 
therein). The unified scheme is further supported by 
the analysis of the radio luminosity functions of the 
postulated parent and aligned populations (Padovani 
& Urry 1992). 

One of the cornerstones of the orientation-based 
paradigm for powerful radio sources has been the 
observation that in the low-frequency 3CRR sam- 
ple (Laing, Riley, & Longair 1983), where the axes 
of the sources should be randomly oriented, the me- 
dian linear extent, £, of the extended radio emission 
from QSRs appears significantly smaller (~50% for 
redshifts z > 0.5 ) than that of RGs (Barthcl 1989). 
However, the lack of such a behavior in the same sam- 
ple at z < 0.5 (Kapahi 1990; Singal 1993a), bolstered 
by similar trends reported recently for a few other 
low-frequency samples a few times deeper in flux den- 
sity than the 3CRR, has evoked serious doubts about 



the unified scheme (Singal 1995, 1996; Kapahi ct al. 

1995) . Here we show that this apparently irrefutable 
inconsistency with the data can be resolved quanti- 
tatively by taking into account the available strong 
independent evidence for temporal evolution in both 
the sizes and luminosities of extragalactic double ra- 
dio sources. 

2. Temporal and Luminosity Dependences of 
Radio Source Properties 

Powerful radio sources are generally accepted to 
have rather limited lifetimes; these are typically es- 
timated to be of the order of 10 7 yr, using the mea- 
sured spatial gradients of spectral index across their 
radio lobes (e.g., Alexander & Leahy 1987; Liu, Poo- 
ley & Riley 1992). Although any particular "spec- 
tral age" determination is uncertain, and probably 
somewhat of an underestimate because of the com- 
plications arising from field inhomogeneities (Siah & 
Wiita 1990; Eilek 1996), inverse Compton losses of 
electron energy to the microwave background radi- 
ation (Wiita & Gopal-Krishna 1990), and backflow 
from the jets into the radio lobes (Liu et al. 1992; 
Schcuer 1995), significant centimeter wave emission 
persisting for times in excess of 10 8 yr is unlikely. Dur- 
ing this active phase a double radio source grows in 
size through a hot circumgalactic medium at a speed 
V which depends weakly on both the radio luminos- 
ity, P {V oc P a , a ~ 0.3) (Alexander & Leahy 1987; 
Liu ct al. 1992) and time (V oc i' 3 , (3 ~ -0.25 to 
0) (Gopal-Krishna & Wiita 1987, 1991; Fanti et al. 

1996) . Further, recent studies (Rcadhead 1995; Fanti 
et al. 1996), combining complete samples of powerful 
radio sources in different ranges of linear size, strongly 
suggest that the growth of a source from a few kpc to 
a few hundred kpc during the source lifetime is accom- 
panied by typically an order-of-magnitude decrease in 
P (largely due to expansion losses), after which the 
nuclear activity subsides rapidly (Rees 1994). Thus, 
we shall take T — 10 7 yr as the characteristic e- 
folding time for the radio luminosity decay, and we 
approximate P oc exp(— t/T), after a short-lived ini- 
tial brightening phase during which the source turns 
on to a level related to the power of its jets. While, 
admittedly, neither fully constrained nor unique, this 
empirically supported characterization of luminosity 
evolution of individual sources appears to give a rea- 
sonable account of the temporal behavior of powerful 
extragalactic radio sources. 
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A key element (Antonucci & Miller 1985; Barthel 
1989; Lawrence 1991) of the unified picture is the crit- 
ical angle ip, which separates the sources classified as 
narrow-line objects (RGs) from those called broad- 
line objects (QSRs). Assuming that all powerful radio 
sources are intrinsically similar, this angle is readily 
determined from the fraction, f q , of sources identi- 
fied as QSRs in a volume-limited sample of sources, 
if the sample is unbiased in orientation. (Alternative 
possibilities, such as the assumption that a subset of 
RGs whose optical spectra show only weak narrow- 
line emission would not exhibit broad lines from any 
viewing direction [e.g., Laing et al. 1994], exist but 
will not be considered further here, as we are inter- 
ested in examining whether the simplest extensions 
of the unification scheme can be reconciled with the 
data.) In practice, such unbiased samples are selected 
at meter wavelengths, where the emission is domi- 
nated by radio lobes presumed to be radiating isotrop- 
ically. Physically, ip is identified with the half angle of 
the polar openings of the dusty obscuring torus sur- 
rounding the active galactic nucleus (e.g., Antonucci 
& Miller 1985; Antonucci 1993; Jaffe et al. 1993). Al- 
though earlier investigations (Kapahi 1990; Lawrence 
1991; Singal 1993a) of the 3CRR sample suggested 
a systematic increase in ip with P (or, z), this re- 
quirement was found to be less compelling in some 
recent studies (Laing et al. 1994; Saikia & Kulkarni 
1994) of the same 3CRR sample. Nonetheless, com- 
parisons of / 9 's found for the 3CRR and three other 
low-frequency complete samples which are selected at 
408 MHz, have an expected z(median) ~ 0.6 (Con- 
don 1993), and arc roughly an order-of-magnitude 
deeper in flux density than the 3CRR sample, dis- 
tinctly indicate a steady increase in tp with P (Gopal- 
Krishna 1995; Singal 1995, 1996). Thus, beginning at 
large values of ~50°-60° for the most luminous radio 
sources (namely, the 3CRR sources at z > 1), ip ap- 
pears to shrink gradually to 20° to 30° for sources in 
the samples selected at 408 MHz which are likely to 
be two orders-of-magnitude less luminous, on average 
(Gopal-Krishna 1995; Singal 1995, 1996). Note that 
we consider the dependence of ip to be primarily on 
(initial) P instead of z, firstly because it is expected 
on theoretical grounds (see below), and secondly be- 
cause such a dependence has also been inferred for op- 
tically selected active galactic nuclei (Falcke, Gopal- 
Krishna, & Biermann 1995, and references therein). 

One likely theoretical explanation for the variation 
in torus opening angle involves the action of ener- 



getic nuclear photons on the dust particles within the 
torus, so that more powerful central engines may be 
expected to have larger values of ip (Nctzer & Laor 
1993; Konigl & Kartje 1994). Note that for a given 
source, ip need not decrease monotonically with time, 
as the secular decay of its radio luminosity (see above) 
can be expected to arise primarily from expansion 
losses (Scheuer 1974; Gopal-Krishna & Wiita 1991), 
rather than from a decline in the intrinsic jet power. 
The diminishing of the radio lobe pressure due to 
the expansion would also lead to a lowering of the 
optical output of the narrow-line emitting filaments 
within and around the lobe, thus explaining the cor- 
relation of the radio luminosity P with the narrow [O 
III] line emission, established by Baum & Heckman 
(1989) and Rawlings & Saunders (1991). The possi- 
bility exists, on the other hand, that ip may increase 
with time, due to a steady rise in the mass of the cen- 
tral engine. The greater gravitational dominance of 
a more massive black hole over the surrounding stel- 
lar core is expected to push outward the inner edge of 
the molecular torus, due to a shift in the region where 
reduced differential rotation allows a piling-up of the 
accreted gas (Yi, Field, & Blackman 1994). Such a 
circumstance would facilitate our explanation for the 
apparent excess of the QSR radio sizes vis-a-vis RGs 
(see below). 

3. The Model 

In light of the above observations and discussion, 
we adopt the following simple relations to broadly 
represent the actual dependences for radio luminos- 
ity, critical angle, and physical (deprojected) size, re- 
spectively: 

P(t) = P exp(-t/T), where T = 10 7 yr; (1) 

iP(P ) = a log^Po/lO^WHz- 1 ), (2) 
where 0.2 radian <ip< ipmax — 1.0 radian; 

L(Po,t)^P*H b2 ; (3) 

where a, b\ and hi are constants (input parameters 
to the model) whose characteristic values can be ap- 
proximated from the empirical estimates for f q , a, 
and P mentioned above. 

We can now compute the RG-to-QSR median lin- 
ear size ratio, R, for a temporally evolving popu- 
lation of sources. Consider a radio luminosity, P, 
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which is below the value for the most powerful radio 
sources created (P ma x)- At any particular time, the 
sources observed at the luminosity level P would in- 
clude not just the young (hence small) sources freshly 
created with that level of luminosity, but also the ag- 
ing sources that were born with higher luminosities in 
the past, but have since then faded down to P, and 
concurrently expanded to larger sizes. Since these 
older expanded radio sources with larger £ would have 
a higher QSR fraction (f q ) arising from their higher 
initial Pq (and correspondingly larger -0), the median 
radio size of QSRs found in intermediate luminosity 
samples with P substantially below P max may well 
approach, or even exceed, that of RGs. This could 
possibly be the explanation for the reported radio size 
"anomaly" mentioned above (Gopal-Krishna 1995), 
which is claimed (Singal 1993a, 1995; Kapahi et al. 
1995) to rule out the orientation based unified scheme. 

To explore this proposal quantitatively, let us as- 
sume that a certain number of sources is "injected" 
continuously at a fixed rate (arbitrarily set to 100 per 
unit time) within a unit volume of space, with radio 
luminosities distributed as 



N(P ) cx P ~ r , for P up to P„ 



(4) 



Our simple model parameterized by Eqns. (1-4) 
can now be used to numerically predict the quasar 
fraction, f q at a given P, as well as the correspond- 
ing ratio, R, of the median projected radio sizes of 
RGs and QSRs. To do this we first need to derive 
the distribution function of ages of sources of a given 
current luminosity P. To do so we start from the full 
distribution function of P and t, which is 



p(P,t)=KP- r exp[-(T-l)t/T}. 



(5) 



^From this we derive the (normalized) distribution 
function of t for a given P to be 

We then generate a large number of sources, the ages 
of which follow this distribution. For each source we 
take its age, U, and its orientation angle from the 
line-of-sight, 0, (randomly selected from to n/2), 
yielding its Po,i (Eqn. 1) and tpi (Eqn. 2). If Q- L < ipi 
then the object is counted as a quasar, otherwise it is 
counted as a galaxy. We then obtain its linear size, 
L (Eqn. 3), and projected linear size, t. By repeating 



this procedure many times, we compute the median 
projected linear size for RGs and for QSRs (yielding 
R), as well as f q , for the adopted value of P. 

In Fig. 1 our simple model is confronted with the 
observations by plotting the combinations of f q and 
R computed as above for different values of P, taking 
some characteristic values for the input parameters 
based on the above discussion and using additional 
observations to set T w 2 in Eqn. 6 (see Dunlop & 
Peacock 1990). In Eqn. 4, P m ax, has been set equal 
to 10 29 W Hz -1 at 1 GHz, characteristic of the most 
luminous radio sources. The runs of the computed 
values of f q and R are shown by three curves, each 
covering a two orders-of-magnitude range in P. Ap- 
proximately the same range is spanned by the data 
points, which are reproduced from Singal (1995), and 
correspond to subsets of the 3CRR, 1-Jy and B3 sam- 
ples. Note that although redshifts for the B3 subset 
are not fully known, its median flux density of ~1 
Jy at 408 MHz suggests (Condon 1993) a z (median) 
^0.6, making it roughly two orders-of-magnitude less 
luminous, on average, than the most luminous 3CRR 
subset (z > 1) shown by the rightmost data point in 
Fig. 1. We further note that (in the absence of com- 
plete rcdshift information) although the size ratios 
given in Singal (1995; also, Fig. 1) actually refer to 
angular sizes, they are not likely to differ significantly 
from the linear size ratios, as most of the sources in- 
volved lie at z > 0.5 where the linear-to-angular con- 
version is relatively insensitive to z (see also Table 
1 in Singal 1993a). Moreover, a similar drop in the 
value of R to around unity near S408 ~1 Jy has been 
inferred (Kapahi et al. 1995) from the linear sizes of 
the Molonglo QSRs and 3CRR RGs. 

As highlighted by Singal (1995, 1996), the data 
points provide a poor match to the prediction of the 
orientation-based unified scheme (shown by the curve 
"U" in Fig. 1) for any single value of ip; even if ip were 
allowed to vary, the predicted trend runs counter to 
the disposition of the data points. It is seen, however, 
that once an allowance is made for the (inevitable) 
temporal evolution broadly characterized by our sim- 
ple, empirically supported model, the predicted f q —R 
profiles for the unified scheme are found to match the 
data points quite well (Fig. 1). 

An interesting implication of the present model is 
that except in the meter wavelength samples selected 
at the highest luminosities, quasars would be system- 
atically older (and physically larger) than RGs. This 
is opposite to the situation envisioned in some alter- 
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native unification schemes (e.g., Hutchings, Price, & 
Gower 1988). 

Finally, a cautionary remark, motivated by our im- 
plicit assumption that broad and narrow-line objects 
can be distinguished at all luminosities and redshifts. 
This assumption may not be valid at lower luminosi- 
ties where any broad component underlying a narrow 
spectral feature may be increasingly hard to detect, 
leading to an underestimate of f q (and, therefore, ip). 
The significance of this potential bias cannot be quan- 
tified using the currently available data, though it is 
unlikely to be a major effect for the relatively lumi- 
nous objects being considered here. 

4. Conclusions 

We emphasize that the temporal evolution incor- 
porated here into the orientational unified scheme is 
exceedingly simple and expressed in terms of just a 
few parameters, all of which are constrained by ob- 
servations. Since this evolutionary scenario is inspired 
by observations and elementary theoretical consid- 
erations, the neglect of this factor in the past at- 
tempts to verify the unified scheme was a major short- 
coming. Consequently, the claimed "mismatch" with 
the radio size data (Singal 1993a, 1995, 1996; Ka- 
pahi et al. 1995) should not mandate dismissal of the 
paramount role assigned to orientation effects in the 
unified scheme. 

The present explanation for the decreasing RG-to- 
QSR size ratio, R, towards lower radio luminosities 
also provides an explanation for an equivalent obser- 
vational result according to which QSRs and RGs ex- 
hibit different l-P correlations. Several authors have 
pointed out that the empirically derived positive l-P 
dependence for RGs (£ oc P x , x ~ +0.3), contrasts 
with practically no l-P dependence found for QSRs 
(x ~ 0), and have argued that this "puzzling" re- 
sult effectively rules out the unified scheme (Chyzy 
& Zieba 1993; Singal 1993a,b, 1995, 1996; Kapahi et 
al. 1995). However, this difference too, can be readily 
understood in our picture. The values of x corre- 
sponding to the three computed model curves shown 
in Fig. 1 arc +0.24 ± 0.07 for RGs and 0.09 ± 0.08 
for QSRs, as deduced from the linear fits to the com- 
puted median sizes at different radio luminosities for 
the three sets of input parameters; for each fit the 
slope x is found to be greater for RGs than for QSRs. 
As a refinement to the present analysis, the use of 
full linear size distributions (rather than the median 



values) would be a useful step. 

In summary, we conclude that the linear size mea- 
surements of radio galaxies and quasars reported up 
to this point do not violate the basic tenet of the 
unified scheme, viz, that radio galaxies are oriented 
closer to the sky plane than are quasars. 

We thank the anonymous referee for suggestions 
which led to improvements in the presentation. This 
work was supported in part by NSF grant AST- 
9102106, NASA grant NAG 5-3098, and Chancellor's 
Initiative Funds at Georgia State University 

REFERENCES 

Alexander, P. & Leahy, J.P. 1987, MNRAS, 225, 1 

Antonucci, R. 1993, ARA&A, 31, 473 

Antonucci, R.R.J. , & Miller, J.S. 1985, ApJ, 297, 621 

Barthcl, P.D. 1989, ApJ, 336, 606 

Baum, S.A., & Hcckman, T. 1989, ApJ, 336, 702 

Chyzy, K.T., & Zieba, S. 1993, A&A, 267, L27 

Cimatti, A., de Serego Alighicri, S., Fosbury, R.A.E., 
Salvati, M., & Taylor, D. 1993, MNRAS, 264, 421 

Condon, J. 1993, in Sky Surveys: Protostarsto Proto- 
galaxies (ASP Conf. Ser. 43), ed. B.T. Soifer (San 
Francisco, ASP), 87 

Draper, P.W., Scarrott, S.M., & Tadhunter, C.N. 
1993, MNRAS, 262, 1029 

Dunlop, J.S., & Peacock, J.A. 1990, MNRAS, 247, 19 

Eilek, J.A. 1996, in Energy Transport in Radio Galax- 
ies & Quasars (ASP Conf. Ser.), eds. P. Hardee, A. 
Bridle, & A. Zensus (San Francisco, ASP), in press 

Falcke, H., Gopal-Krishna, & Biermann, P. 1995, 
A&A, 298, 395 

Fanti, C, Fanti, R., Dallacasa, D., Schilizzi, R.T., 
Spencer, R.E., & Stanghellini, C. 1996, A&A, in 
press 

Gopal-Krishna. 1995, Proc. Nat. Acad. Sci., 92, 
11399 

1996, Proc. IAU Symposium No. 175 (Ex- 

tragalactic Radio Sources), eds. C. Fanti & R.D. 
Ekers (Dordrecht, Kluwer), in press 

Gopal-Krishna & Wiita, P.J. 1987, MNRAS, 226, 531 

1991, ApJ, 373, 325 



5 



Hutchings, J.B., Price, R., & Gower, A.C. 1988, ApJ, 
329, 122 

Jaffc, W., Ford, H., Ferrarese, L, van den Bosch, F., 
& O'Conncll, R.W.O. 1993, Nature, 364, 213 

Kapahi, V.K. 1990, in Parsec-Scale Radio Jets, eds. 
J. A. Zensus, & T.J. Pearson (Cambridge, Cam- 
bridge Univ. Press) 304 

Kapahi, V.K., Athreya, R.M., Subrahmanya, C.R., 
Hunstead, R.W., Baker, J.C., McCarthy, P.J., & 
van Breugel, W. 1995, J. Astrophys. Astr. Supp., 
16, 125 

Konigl, A., & Kartje, J. 1994, ApJ, 434, 446 

Laing, R.A., Jenkins, C.R., Wall, J.V., & Unger, S.W. 
1994 in The Physics of Active Galaxies (ASP Conf. 
Ser. 54), eds. G.V. Bicknell, MA. Dopita, & P.J. 
Quinn (San Francisco, ASP), 201 

Laing, R.A., Riley J.M., & Longair, M.S. 1983, MN- 
RAS, 204, 151 

Lawrence, A. 1991, MNRAS, 252, 586 

Liu, R., Pooley, G., & Riley, J.M. 1992, MNRAS, 257, 
545 

Netzer, H., & Laor, A. 1993, ApJ, 404, L51 

Padovani, P., & Urry, CM. 1992, ApJ, 387, 449 

Rawlings, S., & Saunders, R. 1991, Nature, 349, 138 

Readhead, A.C.S. 1995, Proc. Nat. Acad. Sci., 92, 
11447 

Rees, M.J. 1994, QJRAS, 35, 391 

Saikia, D.J., & Kulkarni, V.K. 1994, MNRAS, 270, 
897 

Scheuer, P.A.G. 1974, MNRAS, 166, 513 

MNRAS, 277, 331 

Siah, M.J., & Wiita, P.J. 1990, ApJ, 363, 411 

Singal, A.K. 1993a, MNRAS, 262, L27 

1993b, MNRAS, 263, 139 

1995, Proc. Nat. Acad. Sci., 92, 11407 

1996, MNRAS, 278, 1069 

Urry, CM., & Padovani, P. 1995, PASP, 107, 803 

Wiita, P.J., & Gopal-Krishna. 1990, ApJ, 353, 476 

Yi, I., Field, G.B., & Blackman, E.G. 1994, ApJ, 432, 
L31 

This 2-column preprint was prepared with the AAS IATfrjX 
macros v3.0. 




0.2 0.4 0.6 0.8 
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Fig. 1. — A plot of R versus f q for the subsets of three 
flux-limited samples selected at meter- wavelengths. 
The data shown with error bars are reproduced from 
Singal (1995). The thin curve "U" depicts the predic- 
tion of the unified scheme for such orientation bias- 
free samples, over a range in the torus half-angle ip. 
The three other curves represent our predictions of 
the unified scheme incorporating temporal evolution, 
as described in the text. Each of the curves covers 
a two orders-of-magnitude range in radio luminos- 
ity, increasing to the right. The input parameters 
(r,a,&i,i/w) are (1.8, 0.38, 0.18, 1.1 rad) for the 
solid curve, (1.8, 0.40, 0.25, 0.9 rad) for the dashed 
curve, and (2.0, 0.40, 0.30, 0.9 rad) for the dotted 
curve. The other input parameters were set (follow- 
ing the empirical results mentioned in the text) to: 
b 2 = 0.8 and -0mm = 0.2 rad. Note that a horizontal 
error bar is also shown for the rightmost data point; 
it represents a typical ler error of these 3CRR data 
points, each of which is based on a relatively small set 
of sources. 
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